Universal behavior of CePdi_ x Rh x Ferromagnet at Quantum Critical Point 
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The heavy-fermion metal CePdi_ x Rh x can be tuned from ferromagnetism at x = to non- 
magnetic state at some critical concentration x c . The non-Fermi liquid behavior (NFL) at x ~ x c is 
recognized by power low dependence of the specific heat C(T) given by the electronic contribution, 
magnetic susceptibility x{T) and volume expansion coefficient a(T) at low temperatures: C/T oc 
x{T) oc a(T)/T oc 1/y/T. We also demonstrate that the behavior of normalized effective mass 
observed in CePdi_ x Rh x at x ~ 0.8 agrees with that of Mjy observed in paramagnetic CeRu2Si2 
and conclude that these alloys exhibit the universal NFL thermodynamic behavior at their quantum 
critical points. We show that the NFL behavior of CePdi_ x Rh x can be accounted for within 
frameworks of quasiparticle picture and fermion condensation quantum phase transition, while this 
alloy exhibits a universal thermodynamic NFL behavior which is independent of the characteristic 
features of the given alloy such as its lattice structure, magnetic ground state, dimension etc. 

PACS numbers: 71.27.+a, 74.20.Fg, 74.25.Jb 



The nature of the non-Fermi liquid (NFL) behavior 
observed in heavy-fermion (HF) metals is still hotly de- 
bated. It is widely believed that the observed behav- 
ior is determined by quantum phase transitions which 
occur at quantum critical points, while the proximity 
of a system to quantum critical points creates its NFL 
behavior brought about by the corresponding thermal 
and quantum fluctuations suppressing quasiparticle exci- 
tations [HHI]. A quantum critical point (QCP) can arise 
by suppressing the transition temperature T c of a ferro- 
magnetic (or antiferromagnetic) phase to zero by tuning 
some control parameters other than temperature, such as 
pressure, magnetic field, or doping x as it takesplace in 
the case of the HF ferromagnet CePdi_ x Rh x or the 
HF metal CeIn3_ x Sn x Q. QCPs are of great interest due 
to their singular ability to influence the thermodynamic 
properties of materials producing the NFL behavior. The 
NFL behavior around QCPs manifests itself in various 
anomalies. One of them is power in T variations of the 
specific heat C(T), thermal expansion a(T), magnetic 
susceptibility x(T) etc. [EUS]. 

Measurements on CePdi_ x Rh x show that around con- 
centration x = x c ~ 0.9 the suppression of the ferro- 
magnetic phase takes place, so that this alloy is tuned 
from ferromagnetism at x = to non-magnetic state at 
QCP with the critical concentration x c - Studies of 
the NFL behavior revealed in the HF metal CePdi_ x Rh x 
0, HI are of great interest since this alloy is a three di- 
mensional ferromagnet. Basing on the theory of critical 
fluctuations which claims that these are responsible for 
the corresponding NFL behavior [J, 0, one can assume 
that the NFL behavior demonstrating by CePdi_ x Rh x 
is to be different from that of CeNi2Ge2 exhibiting a 
paramagnetic ground state [7| or from that of the an- 
tiferromagnetic cubic HF metal CeIn3_ x Sn x [6J. Obvi- 
ously the corresponding critical fluctuations taking place 



at QCPs in the mentioned three different HF metals 
are different, therefore one cannot expect to observe a 
universal behavior demonstrating by these metals, while 
the traditional theory has no grounds to consider these 
QCPs as a single QCP. Moreover, the distinctive features 
between ferromagnetic, antiferromagnetic and paramag- 
netic systems suggest intrinsic differences in their QCPs 
resulting in the difference of their thermodynamic proper- 
ties, and the theory predicts that magnetic, thermal and 
transport properties of these systems have to be different 

At the critical concentration x c , measurements on 
CePdi_ x Rh x show that the specific heat C(T)/T oc 
1/y/T, while around that concentration C/T and x(T) 
coincide in their temperature dependence, C{T)/T oc 
x(T) oc 1/VT [4j, |a]. Moreover, as we shall see it proved 
to be a{T)/T oc 1/VT and this NFL behavior of the 
thermal expansion coefficient coincides with that of a(T) 
observed in the HF metals CeNi2Ge2 |7[ and CeIn3_ x Sn x 
0, HJ. The observed power laws and relationships be- 
tween them can be hardly accounted for within scenarios 
based on the QCP occurrence with quantum and ther- 
mal fluctuations 0, 0, 0, 0] when quasiparticles are sup- 
pressed, for there is no reason to expect that C(T), x(T), 
a{T) and other thermodynamic propreties are affected by 
fluctuations in a correlated fashion. 

These demonstrate that the fluctuations are not re- 
sponsible for the observed behavior, and if they are 
not, what kind of physics determines the NFL behavior? 
Fortunately, the direct observations of quasiparticles in 
CeCoIns have been reported recently |lfj. On the other 
hand, when the electronic system of HF metals under- 
goes the fermion condensation quantum phase transition 
(FCQPT), the fluctuations are strongly suppressed and 
cannot destroy the quasiparticles which survive down to 
lowest temperatures and we can safely suggest that quasi- 
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particles are responsible for the NFL behavior observed 
in HF metals @, [□], H El ■ 

In this letter we show that the NFL behavior of the 
thermal expansion coefficient a{T)/T oc 1/VT observed 
in CePdi_ x Rh x coincides with that of a(T) observed in 
both CeNi 2 Ge 2 and CeIn 3 _ x Sn x . While the NFL be- 
havior of the ferromagnet CePdi_ x Rh x related to the 
uniform temperature dependence of C{T)/T oc x(^) cc 
can be accounted for within the frame- 
work of quasiparticle picture and FCQPT. We demon- 
strate this alloy is of great interest as it exhibits the uni- 
versal NFL thermodynamic behavior at its QCP. This be- 
havior is independent of the characteristic features of the 
given alloy such as its lattice structure, magnetic ground 
state, dimension etc. We also conclude that numerous 
CQPs assumed to be responsible for the NFL behavior 
of the thermal expansion coefficient and other thermody- 
namic properties observed in HF metals can be substi- 
tuted by the only QCP related to FCQPT. 

To study the low temperature universal features of HF 
metals, we use a model of homogeneous HF liquid with 
effective mass M* in order to avoid the complications as- 
sociated with the crystalline anisotropy of solids. This is 
possible since we consider the universal behavior related 
to the power-law divergences of observable values like the 
effective mass, thermal expansion coefficient, specific heat 
etc. These divergences are determined by small (as com- 
pared to those from unit cell of the corresponding recip- 
rocal lattice) momenta transfer so that the contribution 
from larger momenta can be safely ignored. 

To describe the effective mass M* as a function of tem- 
perature and applied magnetic fields B when the system 
approaches FCQPT from the disordered side, x — > x F c, 
we use the Landau equation connecting the effective mass 
M*(T, B) with the bare mass M and Landau interaction 
amplitude F(pi,p 2 ,a;) [14| 
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where n(p, T, B) is the quasiparticle distribution function 

-(e(p,T,B)-»(B)) 



n{p,T,B) = il + exp 



T 



(2) 



Here both the single-particle energy e(p, T, B) and chem- 
ical potential i^(T, B) depend on temperature and mag- 
netic field. It follows from Eq. ([2]) that at B — > and 
T — > the distribution function n(p, T, B) — > 9(p F — p) 
with 9(p F — p ) be in g the step function and we obtain from 



Eq. (HJ) that [3, 

M*(x) = 
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Here N is the density of states of a free electron gas, p F 
is Fermi momentum, F 1 (p F ,p F ) is the p- wave component 



of Landau interaction amplitude, A and B are constants. 
Since Landau Fermi liquid (LFL) theory implies the num- 
ber density in the form x = p F /3ir 2 , we can rewrite the 
amplitude as F 1 (p F ,p F , %) = F 1 (x). When x — > x F c, 
F 1 (x) being a function of x achieves some value at which 
the denominator tends to zero so that the effective mass 
diverges at T = as seen from Eq. ([3]). 

At first let us consider the dependence of the effective 
mass on temperature. Upon using Eq. ([3]) and intro- 
ducing the function Sn(p, T) — n(p,T) — 9{p F — p), we 
transform Eq. (fTJ and it takes the form 
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(4) 

We integrate the second term on the right hand side of 
Eq. |4]) over the angle variable and use the notation 



Fi (p F ,p, x) = Mp F J p F 



8F(pf,P,x) dSl 



dpi 



(2tt) 3 ' 



(5) 



and substitute the variable p by z, z — (s(p) — /i)/T . 
Since in HF metals the band is flat and narrow, we use 
the approximation (e(p) — /i) ~ Pf{p — Pf)/M*(T) and 
taking into account Eqs. ((4]) and ([5|) finally obtain 
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Here the factor a.\ = TM*(T)/p F . The Fermi momen- 
tum p F is defined from the relation e{p F ) — /i. We first 
assume that M*(x) is finite and a\ <^ 1. Then upon 
omitting terms of the order of exp(— 1/ai), we expand the 
upper limit of the second integral on the right hand side 
of Eq. (|6|) to oo and observe that the sum of the second 
and third terms represents an even function of a\ . These 
are the typical expressions with Fermi-Dirac functions as 
integrands and can be calculated using standard proce- 
dures [13] . We conclude that at T < T F - p%/M*(x) the 
sum represents a T 2 -correction to M*(x) and the system 
demonstrates the LFL behavior [l8| . When x — > x pc the 
effective mass diverges and both T F — > and l/o^ — > 0, 
while the temperature interval over which the LFL be- 
havior takes place is vanishing In that case, Eq. 
([6]) becomes homogeneous and the second integral on the 
right hand side can be omitted. As a result, we can esti- 
mate that 



M*(T) oc 



(7) 



Equation ^ shows the universal power low behavior of 
the effective mass which does not depend on the inter- 
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particle interaction. We illustrate this behavior by calcu- 
lations using a model functional 

E[n(p)] = [ — — 7"—tt + - / V(pi - p 2 ) 
L yPU J 2M (2tt) 3 2 J VP1 V2 ' 



(27T) 

x n(pi)n(p 2 ) 



rfpirfp 2 

(2tt) 6 ' 

with the inter-particle interaction 

cxp(-/3 |p|) 



V(p) = <?o- 



(8) 



(9) 



We normalized the effective mass by M, M* — 
M*(T)/M, temperature T by the Fermi energy e%, 
T = T/e'p and use the dimcnsionless coupling constant 
g = {g M)/(2ir 2 ) and f3 = (3 a p F . FCQPT takes place 
when the parameters reach their critical values, = b c 
and g = g c , in our case b c = 3 and g c = 6.7176. In FigQ] 



On the other hand, the specific heat 
C(T) 



dT 



(12) 



As a result, at T — > the Griineisen ratio T(T) tends 
to some constant value rather than diverges as in the 
case when the electronic system is on the ordered side of 
FCQPT dUEll 



r(T) 



<*(T) 
C{T) 



— const. 



(13) 



Since the magnetic susceptibility x(T) cx M*(T) and 
both the Sommerfeld coefficient C(T)/T cx M*{T) and 
a(T)/T cx M*{T) we conclude that at T -> 

^ cx X (T) cx ^ cx M*(T) cx JL. (14) 
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FIG. 1: Calculated entropy S(T) as a function of tempera- 
ture. Continuous line represents S(T) — ciyT, where c\ is a 
parameter. Solid squares are the results of calculations based 
on functional ((8]). 

the evolution of the low temperature entropy is shown. 
The calculated behavior of S(T)/T cx M*(T) cx l/VT is 
in accord with Eq. |7]). 

Now consider the thermal expansion coefficient a(T) 
given by [17 1 



a(T) 



1 / a(lnV) 
3 dT 



1 W/s) 
3V I 



(10) 



Here, P is the pressure and V is the volume. The com- 
pressibility K{x) is not expected to be singular at FC- 
QPT and is approximately constant Inserting into 
Eq. ([ID} the entropy S(T) cx VT, we find that 



a(T) 



M*T 



T. 



(11) 



FIG. 2: The thermal expansion coefficient a(T) as a function 
of temperature in the interval 100 mK < T < 6 K. Continuous 
curves are fits for x = 0.87 and x = 0.90 data [J] based on Eq. 
([7]) and represented by function a(T) = civT with ci being 
a fitting parameter. 

At this point, we consider how Eq. and the be- 

havior of the effective mass given by Eq. ([7]) correspond 
to experimental observations obtained on CePdi_ x Rh x . 
Measurements of the thermal expansion coefficient a(T) 
on CePdi_ x Rh x with x = 0.87 and x = 0.90 0] are 
shown in Fig. [2l It is seen that the approximation 
a(T) — c\\pf for composition x — 0.90 is in good agree- 
ment with facts over two orders of magnitude in the tem- 
perature range from 6 K down to 100 mK, and measure- 
ments on CeNi 2 Ge2 and CeIn 3 _ x Sn x [|| demonstrate 
the same behavior. While CePdi_ x Rh x is a three di- 
mensional ferromagnet 0,3, CeNi 2 Ge2 exhibits a para- 
magnetic ground state [jj and CeIn3_ x Sn x is antiferro- 
magnetic cubic metal @ . We conclude that the observed 
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uniform behavior of the thermal expansion coefficient of 
these metals is determined by quasiparticles and FCQPT 
rather than by different magnetic quantum critical points 
and corresponding fluctuations. Measurements of the 
specific heat C(T) on CePdi_ x Rh x with x — 0.87 and 
x = 0.90 show a power law T dependence. These are de- 
scribed by a C(T)/T = AT~ q formula with the exponent 
q ~ 0.5 — 0.4 and A is a constant, around that concentra- 
tion C(T)/T and x{T) coincide in their T~ q temperature 
dependence 0, 0] • We conclude that the results given by 
Eq. (fT4")) agree with facts. 
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FIG. 3: Normalized magnetic susceptibility 

Xac(T,B)/xac(Tm,B) for CeRu2Si2 in magnetic fields 
0.20 mT (squares), 0.39 mT (upright triangles) and 0.94 
mT (circles) is plotted against normalized temperature 
T/Tm 24!\. The susceptibility reaches its maximum value 
Xac(Tm , B) at T = Tm. Normalized 4f electron contri- 
bution (AC(T)/T)/(AC(Tm)/T m ) to the specific heat of 
CePdi_ x Rh x with x = 0.80 versus normalized temperature 
T/Tm is shown by downright triangles Here Tm is the 
temperature at the peak of AC(T)/T. The solid curve 
traces the universal behavior of the normalized effective mass 
determined by Eq. (I18|l . 

Consider the case when the concentration x deviates 
from the critical value xfc so that the system is moved 
to the disordered side of FCQPT and temperature Tp be- 
comes finite. As a result, at T <C Tp the system exhibits 
the LFL behavior with the effective mass being approx- 
imately constant, M*(T) ~ M*(x) + ai(T/T F ) 2 , where 
the term a\(T/Tp) 2 represents correction to M*(x). If 
x — xfc then the application of magnetic field making 
Tp finite restores the LFL behavior with the effective 



mass depending on B as [22|, |23 1 

M*(B)cx(B-B c0 y 



-2/3 



(15) 



where B c q is the critical magnetic field which drives both 
a HF metal to its magnetic field tuned QCP and the 



corresponding Neel temperature toward T — 0. In some 
cases B c q = 0, for example, the HF metal CeRu2Si2 is 
characterized by B c0 — and shows neither evidence of 
the magnetic ordering, superconductivity nor the LFL 
behavior down to the lowest temperatures [24J . 

At T ~ Tpj the effective mass depends mainly on tem- 



perature 



M*(T) cx T- 2/3 . 



(16) 



Then, at elevated temperatures T Tp, the behavior 
of the effective mass is given by Eq. © [18]. Therefore 
at T < Tp the behavior of the effective mass can be 
described by a simple function 25] 



1 



c 2 y 



M*(B,T) 
M*{B) ~ 1 + c 3 ?/ 8 / 3 : 



(17) 



which represents an approximation to solutions of Eq. 
|T|) that agrees with Eqs. (fl~5")) and (fTrJ]) . Here y = 
(T I (B — B c o)), C2 and C3 are fitting parameters. Since the 
effective mass reaches its maximum value at some 
y = yM we define a normalized effective mass 

as M* N (T, B) = M*{T, B)jM* M . Taking into account Eq. 



(fT7|) and introducing the variable z 
the function 



M* N (z) 



1 1 + c 2 z 2 



M* M 1 



c 3 z 



y/yM we obtain 



(18) 



8/3 ' 



which describes a universal behavior of the effective mass 
M* N (z). In the case of finite M*(x), Eq. (B) is valid 
at T ~ Tp if M*(T,B)/M*(x) < 1 because the term 
1 /M* (x) on the right hand side of Eq. (01 being a small 
correction to the effective mass can be omitted (2(|. It 
is seen from Eq. (|18l) that M%(z) reaches its maximum 
value at z = 1, M^(z = 1) = 1. 

The effective mass M*(T,B) can be measured in ex- 
periments on HF metals. For example, as it follows from 
Eq. CL1 M*(T,B) cx C/T and M*(T,B) cx xac(T) 
where xac(T) is the magnetic susceptibility. If the cor- 
responding measurements are carried out at fixed value 
of magnetic field B (or at fixed value of the concentra- 
tion x and B = 0) then as it follows from Eq. (fT7|) the 
effective mass reaches the maximum at some tempera- 
ture Tm- Upon normalizing both the effective mass by 
its peak height at each field B and the temperature by 
Tm, we observe that all the curves should demonstrate 
a scaling and collapse on the single curve given by Eq. 

As shown in Fig. [3l the behavior of the normalized 
susceptibility X ac( z ) = Xac(T/T m ,B)/ X ac(1,B) = 
Mpf(z) obtained in measurements on the HF paramag- 
netic CeRu2Si2 [24j is in accord with the approximation 
given by Eq. (jTHJ) - Since the crossover temperature T* 
from the regime given by Eq. (|16p to the regime given 
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by Eg ([7]) is proportional to the magnetic field, T* oc B 
3 2M i we expect that the temperature range over which 
the scaling takes place shrinks when the applied magnetic 
field B is diminished. It is seen from Fig. [3] that the de- 
viation of the data corresponding to the smallest value 
of the magnetic field B = 0.20 mT (shown by squares) 
is largest at the elevated normalized temperature, while 
the slope of this curve tends to that of curve described by 
Eq. |7|). At small temperatures as seen from Fig. [3J both 
the curve given by Eq. (|18|) and the effective mass de- 
termined by Eq. (TTS"]) agree perfectly with facts collected 
in measurements on CeRu2Si2 whose electronic system 
is placed at FCQPT [3]. As to the normalized 4f con- 
tribution (AC(T)/T)/(AC(T M )/T M ) = Mn( t ) (shown 
by downright triangles in Fig. [3]) to the specific heat of 
CePdi_ x Rh x with x = 0.80 [5j, the scaling takes place up 
to relatively large temperatures because the deflection of 
the x — 0.8 from the critical concentration xfc — 0.9 is 
big, while T* cx B oc \xfc — x\ [26]. As a result, at di- 
minishing temperatures the scaling is ceased at relatively 
high temperatures as soon as the LFL behavior sets in. 

In summary, we have shown that the NFL behavior 
of the thermal expansion coefficient a(T) observed in 
CePdi_ x Rh x at the critical concentration x c ~ 0.9 co- 
incides with that of a(T) observed in both CeNi2Ge2 
exhibiting the paramagnetic ground state and antifer- 
romagnetic cubic HF metal CeIxi3_ x Sn x . We have also 
shown that the behavior of the normalized effective mass 
observed in CePdi_ x Rh x at x ~ 0.8 agrees with that 
of observed under the application of magnetic field 
in paramagnetic CeRu2Si2 and concluded that these al- 
loys exhibit the universal NFL thermodynamic behavior 
at its QCP. The outlined behavior is independent of the 
characteristic features of the given alloys while numerous 
CQPs assumed to be responsible for the NFL behavior of 
different HF metals can be substituted by the only QCP 
related to FCQPT. 

This work was supported in part by RFBR, project No. 
05-02-16085. 
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